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Consider a multichannel closed ring with disorder. In the semiclassical treatment its conduc- 
tance is given by the Drude formula. Quantum mechanics challenge this result both in the limit 
of strong disorder (eigenstates are not quantum-ergodic in real space) and in the limit of weak 
disorder (eigenstates are not quantum-ergodic in momentum space). Consequently the analysis of 
conductance requires going beyond linear response theory, leading to a resistor network picture of 
transitions between energy levels. We demonstrate that our semi-linear response theory provides a 
firm unified framework from which the "hopping" phenomenology of Mott can be derived. 



I. INTRODUCTION 

The theory for the conductance of closed mesoscopic 
rings has attracted a lot of interestii^i^iii^i^iaiSiiSiiiii^. 

In a typical experiments^ a collection of mesoscopic rings 
are driven by a time dependent magnetic flux <I>(t) which 
creates an electro-niotive-force (EMF) — $ in each ring. 
Assuming that Ohm's law apphes, the induced current is 
/ = — G<f> and consequently Joule's law gives 

yV = G<i>^ — rate of energy absorption (1) 

where G in this context is called the conductance. For 
diffusive rings the Kubo formula leads to the Drude for- 
mula for G. A major challenge in past studies was to cal- 
culate the weak localization corrections to the Drude re- 
sult, taking into account the level statistics and the type 
of occupationiSiiiiii^. It should be clear that these cor- 
rections do not challenge the leading order Kubo-Drude 
result. 

One wonders what happens to the Drude result if the 
disorder becomes weak (ballistic case) or strong (Ander- 
son localization case) . In both cases the individual eigen- 
function become non ergodic: a typical eigenfunction do 
not fill the whole accessible phase space. In the ballistic 
case a typical eigenfunction is not ergodic over the open 
modes in momentum space, while in the strong localiza- 
tion case it is not ergodic over the ring in real space. 

Lack of quantum ergodicity implies that the pertur- 
bation matrix is very structured and/or sparse. Conse- 
quently the calculation of G requires a non-trivial exten- 
sion of linear response theory (LRT) . Such extension has 
been proposed in Ref.^** and later termed "semi linear 
response theory" (SLRT)!^. Within SLRT it is still as- 
sumed that the transitions between levels are given by 
Fermi-golden-rule, but a resistor network analogy is used 
in order to calculate the overall absorption. In order 
to have non-zero absorption we must have connected se- 
quences of transitions. Thus the calculation of the energy 
absorption in Eq. ([l} is somewhat similar to solving a per- 
colation problem. The "percolation" is in energy space 
rather than in real space. 

In a previous work we have worked out results for G in 
the case of a ballistic ring^^. In this work we would like 
to explore the other extreme case of strong disorder. If 
we go to the literature we find a strange twist. For the 



AC conductance Moti^il (see also^s) has used the Kubo 
formula in order to predict ~ w^j log(aj)|''+^ dependence 
of the conductivity, where d — 1,2,3 is the dimension- 
ality of the sample. On the other hand for the DC cal- 
culation Mott has abandoned the Kubo formalism and 
has adopted a phenomenological variable range hopping 
(VRH) picturei^, which can be regarded as an approx- 
imation for a more elaborated (but still phenomenolog- 
ical) resister network picture^SiSiiS^. The "ad hoc" ap- 
proach to hopping is obviously bothering. For example it 
is not clear whether the effect of low frequency noisy driv- 
ing should treated like "low temperature" or like "small 
frequency"—. The main purpose of this Letter is to ex- 
plain that the hopping picture is a natural outcome of 
SLRT, and hence can be regarded as a natural exten- 
sion of LRT. This automatically resolves such conceptual 
problems, and opens the way for further refinements of 
the theory. 

The outline of this Letter is as follows: We summarize 
the LRT and the SLRT recipes for the calculation of the 
conductance. Then we demonstrate that all the known 
results for the conductance can be derived from the same 
theoretical framework, without any extra assumptions. 



II. THE MODEL PARAMETERS: 

We consider a ring of length L and mean free path 
^. The particles are non-interacting "spinless" electrons 
with charge e. The Fermi energy and the Fermi velocity 
are Ep and Wp respectively. The one-particle density of 
states at the Fermi energy is 
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The number of open modes A4 reflects the cross section 
of the ring. The GeometricFactor depends on the dimen- 
sionality d = 1,2,3 and equals unity for d — 1 network 
systems. In what follows, for sake of presentation, and 
without loss of generality, we set in all expressions the 
d = 1 geometric factor, and use units such that h ~ 1. 

The mean level spacing is A = g~^. This is the "small" 
energy scale. We also have the "Thouless energy" or 
its equivalent which we denote as A;,. For ballistic ring 
Afe = A^A is associated with the time scale L/vp, while 
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in the diffusive regime Af, = (£/L)A1A is associated with 
the ergodic time. In the strong disorder regime we shaU 
see that the relevant time scale for our analysis is the 
breaktime t* which is related to the localization length £^ . 

Within the framework of both LRT and SLRT it is es- 
sential to realize that the combined effect of the driving 
and the environment is to "broaden" the energy levels 
and to induce "relaxation". Accordingly we distinguish 
between F which can be interpreted as the dephasing or 
decoherence rate (analogous to I/T2 in NMR studies), 
and 7,ix which is the relaxation rate (analogous to 1/Ti 
in NMR studies). We assume that F is much larger than 
A but much smaller compared with Af,. We also assume 
that the driving source is "essentially DC driving". This 
means that the driving frequency is much smaller com- 
pared with A;,. LRT applies whenever the driven transi- 
tions are much slower compared with the relaxation rate. 
SLRT is an extension of LRT that applies in the oppo- 
site circumstances, namely if the relaxation process can 
be neglected during the time that the absorption rate get 
stabilized. This is explained in length ini^, and will be 
emphasized again in a later paragraph. 

III. THE SEMICLASSICAL LRT 
CALCULATION: 

The semiclassical Kubo formula once expressed with 
the velocity-velocity correlation function is 

G = ^^. x-(-) J {v{t)v{0))dt (3) 

where v is the velocity of the particle along the ring. For 
the purpose of later reference we note that Eq. ^ [with 
Eq. ([3]) substituted] is easily generalized to the case where 
the driving is noisy. Instead of a multiplication by $^ we 
have in general an integration over the all the Fourier 
components of the driving. Namely, 

^ = (^-^2il) J ^'^^H^H (4) 

where C'(lu) is the Fourier transform of the velocity- 
velocity correlation function, and F{uj) is the power spec- 
trum of $(t). 

The Kubo formula parallels the expression for the spa- 
tial diffusion coefficient, which is the integral over the 
velocity-velocity correlation function: 

'^ = 11^ HtHO))dt (5) 

In the case of pure DC driving {F{uj) = ^^S{uj)) one ob- 
tains the Einstein relation: 

G = g.[ffv (6) 

If the driving is noisy, then in general G > gp{e/ L)^T>. 
This reflects the simple observation that an insulator 
{V ~ 0) may have a finite AC conductance (finite W) . 



IV. THE QUANTUM LRT CALCULATION: 

Eq.Q is formally valid also in the quantum mechan- 
ical case, and can be regarded as the outcome of the 
Fermi-Golden-rule picture. However one should use the 
proper expression for C{uj), taking into account the levels 
statistics and the level broadening: 

CU^) - R{oj)C^,{cu) (7) 
Ciu) « SriLo)*CU^) (8) 

where i?(w) takes care for the level statistics, ^r(i^) is the 
line shape of the level broadening, and * indicates a con- 
volution. It should be clear that with pure DC driving 
Eq.(|4]) gives zero due to the discreteness of the energy 
spectrum, unless there is some mechanism that "broad- 
ens the levels" so as to have effectively a contimmm. It is 
customary to express the quantum version of the Kubo 
formula using the matrix elements of the velocity opera- 
tor: 

G = Trh(^^y |i-„„P SriEn-Ep) 5v{E„-E„) (9) 

n^ra 

The T broadened delta should be interpreted as the 
derivative of the Fermi occupation function, while the 
functional shape of the F broadened delta function will 
be discussed later. 

V. THE QUANTUM SLRT CALCULATION: 

A loose way to write the last expression is 

G = T,h[^)\l{{\v^,r^?)) (10) 

The interpretation of ((...)) within the framework of the 
tradition LRT is implied by Eq. ^ . It involves the spec- 
ification of the level broadening parameter F and of the 
occupation temperature T. 

The implicit assumption in LRT is that the driving 
induced transitions are much slower compared with the 
environmentally induced relaxation. In other words it is 
assumed that the relaxation is very effective in "killing" 
any quantum effect that goes beyond first order perturba- 
tion theory. SLRT, unlike LRT is aimed in taking account 
also the opposite circumstances in which the possibility 
to make connected sequences of transitions becomes es- 
sential in order to get absorptioni^. FoUowingi^ii^ we 
regard the energy levels as the nodes of a resistor net- 
work. We define 

g,„„ = 2g-^ 5r{E„,~Er.) (11) 

Then it is argued that {{\vnm\'^))^^ is the resistivity of 
the network. It is a simple exercise to verify that if all 
the matrix elements are the same, say |wnmP = cr^, then 
((knmP)) = too. But if the matrix is structured or 
sparse then ((|wnmP)) is in general much smaller com- 
pared with the RMS value of the matrix elements. 
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VI. LEVEL BROADENING: 



VII. MODERATE DISORDER: 



Neither LRT nor SLRT are self contained without the 
specification of the level broadening. For the purpose 
of analysis we can regard the environmental fluctuations 
as a noisy driving source. This driving source induces 
decoherence: It can be argued (see for example Refj^) 
that the effect of decoherence is to multiply the velocity- 
velocity correlation function by an exponential factor 
exp(— r|t|). This means that 



(5r(t^) 



1 



TT w2 + r2 



(12) 



It can be further argued^i that an estimate of the rate T 
can be obtained using the Fermi golden rule (FGR): 



r 



(13) 



where F{uj) is the power spectrum of the environmentally 
induced noise. At low temperatures (see later) the noise 
power spectrum is effectively very narrow. Schematically 
we can write F(uj) — e^STi^^), where e is the noise am- 
plitude, and T is the temperature. 

The FGR expression Eq. lfTSl) parallels Eq.([4]), and can 
be expressed using the matrix elements of the perturba- 
tion as in Eq. ([9]) , or it can be written loosely in the style 
of Eq.Hni): 



r = £' 



X27rh(^^y {{\Vnrrf)) (M) 



The latter version is common in elementary textbooks 
(the FGR calculated decay rate is proportional to the 
squared matrix element times the density of states). It 
is implicit in this expression that the FGR transitions 
5{Em—En) have a width (broadening) which is deter- 
mined by the normalized power spectrum of the environ- 
mental fluctuations. 

So far we simply re-stated the standard textbook 
version of FGR using our notations. But it should 
be clear that the naive application of the FGR recipe 
is as problematic as the the naive application of the 
Kubo formula^^: the first order (transient) transitions 
to "neighboring" states are not enough in order to get 
a non-transient decay. Rather, the possibility to have 
connected sequences of transitions is essential in order to 
have a long time decay. Thus we conclude that an analo- 
gous SLRT recipe, with gnm defined in a similar fashion 
as in Eg. ljlip . should be applied in order to determine 
the long time value of F. We point out again that in the 
calculation of F the role of the broadened delta is played 
by the normalized power spectrum of the environmental 
fiuctuations. This means that at low temperatures (see 
further discussion in a later paragraph) we have in the 
F oriented version of Eq. ITTI) a thermally broadened delta 
function 5T{Em—En). 



The effect of hard chaos or disorder is to randomize the 
velocity. Within the framework of the semiclassical LRT 
calculation we always get the Drude formula irrespective 
of the strength of the disorder. Following Drude it is 
customary to write 



{v{t)viO)) 



vl exp 



\t\ 



or equivalently 



This leads to the Drude result 



2'Kh L 



(15) 



(16) 



(17) 



Once we turn to the quantum calculation, we have to be 
more careful. For moderate disorder the eigenfunctions 
are ergodic, and therefore the distinction between the 
LRT recipe and the SLRT recipe is not important. Fur- 
thermore, using a random wave conjecture for the eigen- 
states, one recovers the Drude expression C^i[lo) as an 
approximation for the quantum C{oj). In order to do 
a better job, the spectral function R{Ld) is introduced. 
This function takes into account the level statistics: for 
large uj it equals unity, while for small oj it refiects the 
repulsion between levels. The "level broadening" effect 
is taken care of by the convolution with (Sr(^), as in- 
dicated in Eq.®. The standard LRT regime is hav- 
ing A ^ F ^ Afc. In this regime the introduction of 
R{uj) implies so-called weak localization corrections to 
the Drude result. These are found to be of order A/F. 
See RefsiiSiiiii^, and alsc^^ for the "quantum chaos" 
point of view. 



VIII. WEAK DISORDER: 

The SLRT recipe for the calculation of G can be re- 
garded as an extension of the LRT recipe. The results of 
SLRT become very different from those of LRT once the 
perturbation matrix Vnm is either structured or sparse. 
The case of weak disorder (ballistic case) has been ana- 
lyzed in a previous Letter—. In the non trivial ballistic 
regime (1 ^ t/L ^ M.) each eigenfunctions occupy a 
large but finite fraction of open modes. One may say 
that the eigenfunctions are "localized" in mode space. 
This lack of quantum-ergodicity implies sturctures and 
sparsity. Gonsequently the SLRT result is not merely a 
small weak localization correction: the leading order re- 
sult is no longer Drude. In what follows we address the 
other extreme case of strong disorder (Anderson localiza- 
tion case). Also here we are going to see that the leading 
order result is not Drude. 
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IX. STRONG DISORDER: 

The first step is to figure out how R{uj) look like for 
a system with strong localization. The initial spread 
of a wavepacket is diffusive with {{x{t) — x(0))^) = 2I?o^ 
where Vq = v^£. But for long time {{x{t) - x{0)f) = ij, 
where £^ is the localization length. This implies a break- 
time at t* =£|/X'o: and therefore the velocity- velocity 
correlation function is modified in the frequency range 
< (l/t*). Namely, the above implies 



1 



1 + {t*Lu)- 



(18) 



Further considerations which are based on Mott's picture 
of resonances imply an extra |loga;|''+^ factor at small 
frequencies. We can argue in advance that this type of 
correction has no importance in SLRT because it reflects 
a very sparse contribution to the perturbation matrix, 
that cannot lead to connected sequences of transitions. 

Even if we eliminate Mott's resonances, still the low 
frequency behavior of C{u!) reflects a very sparse ma- 
trix |w„mP- Within the framework of SLRT the sparse 
component of |w„rnP does not contribute to the diffu- 
sion. Only the non-sparse component allows connected 
sequences of transitions. We argue that the non sparse 
component is obtained by multiplying Cci(w) by 



exp 



l/d' 



(19) 



The reasoning is as follows: The matrix elements of 
states n that are located within range \xn ~ x„i \ < r from 
a given state m constitute a connected grid with spac- 
ing Ar = [L/rY/S. and typical value \vnm\ oc exp(— r/^^). 
From the equation Ar — lu we deduce that the "volume" 
of states that contribute a connected grid for lo transi- 
tions has a radius r = {A/loY^'^L. Hence we deduce the 
above formula, where = {L/£^YA. 



X. THE HOPPING PICTURE: 

Assuming that the coherence time ~ l/F is much 
smaller compared with the breaktime t* we observe that 
the du! integral of Eq.(|3]) with (I7][5]) and the Lorentzian 
of (fT^ is dominated by the tail [uj > (l/t*)] leading to 
the result 



V 



Tt*VQ 



(20) 



This is as expected from heuristic considerations. It de- 
scribes a random walk hopping process with steps of size 
and time x^. 

The issue is to obtain an explicit expression 
for F. Thermal noise is "white" at high temperatures 
{F{lo) oc T), with a very large temperature-independent 
cutoff frequency. In such case Ea. p3p implies F cx T 



which is not very interesting. Low temperature noise, un- 
like white noise, has an exponentially decaying emission 
tail, and consequently, it has effectively a very narrow 
span of frequencies: 



exp 



H 
T 



(21) 



In the above expression we have neglected term whose 
exponent depends of the detailed spectral properties of 
the bath. The calculation of F with Eg.pS)) involves a 
duj integral over 



exp 



T 



X exp 



l/d 



(22) 



This is mathematically equivalent to the VRH integrali^. 
There the optimization is over the range of the hopping r, 
while here it is over the associated frequency lu. The op- 
timal frequency is cx T'^/^'^^^^ leading to the VRH es- 
timate V oc exp[(-To/T)i/(i+'*)] where Tq is a constant. 
It should be clear that the VRH estimate is an approxi- 
mation for the resistor network calculation with Eg. pTj) . 
The VRH estimate works quite well for d > 1, but gives 
the wrong exponent (T~^/^) for d — 1. The correct ex- 
ponent (T^^) in the latter case is implied by the absence 
of percolation. 



XI. CONCLUSIONS: 

We have established that SLRT provides a firm unified 
framework for the calculation of the conductance. Ad 
hock phenomenology is not required in order to estab- 
lish the resistor network "hopping" picture, from which 
Mott's VRH approximation is derived. It should be clear 
that the generalized resistor network picture of SLRT is 
not limited to the strong disorder regime: it allows on 
equal footing the calculation of the conductance in the 
other extreme case of ballistic motion. The importance 
of this approach is also in its potential capabilities: being 
a natural extension of LRT it also allows, in principle, the 
incorporation of many-body effects. In the latter case the 
calculation of conductance is reduced, as in LRT, to the 
analysis of the matrix elements of the current operator, 
whatever are the interactions involved. 

The long standing puzzle regarding "AC conductance" 
versus "DC conductance" that has been discussed in the 
introduction is automatically resolved by our theoretical 
framework. We claim that there is no crossover from 
"AC conductance" to "DC conductance" as a function 
of frequency. Rather, we argue that there is a crossover 
from "AC conductance" to "DC conductance" as a func- 
tion of the driving intensity. Using the terminology of 
Ref.~, the crossover is from the "spectroscopic" to the 
"mesoscopic" result for the conductance: The Kubo for- 
mula of LRT applies to a spectroscopic measurement of 
the conductance, where the driving is assumed to be very 
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weak compared with the relaxation processes; The hop- 
ping picture that emerges from SLRT appUes in the meso- 
scopic regime where the relaxation is assumed to be slow 
compared with the driving-induced transitions. 
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